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Abstract 

We present a new code, RCF(”Radiative-Collisional code based on FAC”), which 
is used to simulate steady-state plasmas under non local thermodynamic equilibrium 
condition, especially photoinization dominated plasmas. RCF takes almost all of the 
radiative and collisional atomic processes into rate equation to interpret the plas¬ 
mas systematically. The Flexible Atomic Code (FAC) supplies all the atomic data 
RCF needed, which insures calculating completeness and consistency of atomic data. 
With four input parameters relating to the radiation source and target plasma, RCF 
calculates the population of levels and charge states, as well as potentially emission 
spectrum. In preliminary application, RCF successfully reproduces the results of a 
photoionization experiment with reliable atomic data. The effects of the most impor¬ 
tant atomic processes on the charge state distribution are also discussed. 
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1. Introduction 

Non local thermodynamic equilibrium (NLTE) exists in a wide variety of astrophysical 
and laboratory plasmas. Examples of NLTE astronomical plasmas are the stellar corona, inter¬ 
stellar nebulae and some other low density ionized plasmas. X-ray satellites, such as Chandra 
and XMM-Newton, provided large amount of high resolution spectra from astronomical objects, 
many of which are in NLTE. In laboratory, NLTE exists in laser produced plasmas, tokamaks 
and Z-pinch based experiments. 

In the present paper, we introduce a new NLTE plasma computer code, which we called 
Radiative-Collisional code based on FAC (abbreviated RCF). RCF is a radiative-collisional 


1 


Table 1. Atomic processes in RCF (Salzmann 1998). 


Reaction 

Direct Process 

Inverse Process 

Xi jf “h hiy 

Spontaneous Decay 

Photo excitation 

Xi,j -t e ^ W j' + e' 

Electron Impact Excitation 

Electron Impact Deexcitation 

X{ j hv X— j' “h c 

Photoionization 

Radiative Recombination 

+ e" 

Electron Impact Ionization 
Autoionization 

Three-Body Recombination 
Dieletronic Captrue 


code that includes photoionization as well, thereby it is appropriate especially to astrophysical 
plasmas. In the following we show its accuracy relative to other codes. 

Several similar codes, which are being used for analysis of astrophysical spectra, have 
already been published in the literature. Examples are GALAXY (Rose 1998, Foord et ah 
2004, Foord et ah 2006, Rose et ah 2004), NIMP (Djaoui & Rose 1992, Rose et ah 2004), 
FLYCHK (Chung et ah 2003), CLOUDY (Ferland et ah 1998), XSTAR (Kallman et ah 1996, 
Kallman et ah 2004, Kallman & Bautista 2001, Bautista & Kallman 2001, Boroson et ah 2003), 
PhiCRE (Salzmann et ah 2011, Wang et ah 2011) and SASAL (Liang et ah 2014). Some of 
these codes are used to interpret laboratory experiments (GALAXY, FLYCHK, NIMP and 
PhiCRE), while others are used in analysis of astrophysical spectra (CLOUDY and XSTAR), 
while RCF is designed to be applicable to both of the conditions above. 

The aim of this paper is to give a detail introduction of RCF and to present its application 
to a phoionization experiment. Section 2 presents the model and the method of to calculate 
the atomic rate coefficients. In section 3, we apply RCF to reproduce the iron charge state 
distribution of a photonionization experiment, together with a discussion of the importance of 
the various atomic processes. A short summary is given in the last section. 

2. Rate equation and Atomic Data 

2.1. Rate Equation 

RCF is a steady-state collisional-radiative optically-thin model. Its rate-equation 
(Salzmann 1998) is 

dNi j 

——^ = populating processes — depopulating processes = 0 
dt 

where Nij is the density of the jth level of the Ah charge state. The processes included are 
the ionization and recombination between neighboring charge states and excitation and de¬ 
excitation within the same charge state. Their inverse processes are also taken into account by 
detailed balance principle. The processes included in Eq.(l) are listed in Table 1. 
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2.2. Atomic Data and Reaction Rates 


The atomic data for RCF are calculated by FAC (Gu 2008). FAC is a fully relativistic 
software package that computes various atomic data, which has been widely used in astrophys- 
ical and laboratory research (Gu 2008). With the ions and their configurations as input, the 
SFAC interface can supply j — j coupled energy levels (Eij), bound-bound spontaneous decay 
rates bound-bound electron-collision excitation (CE) cross sections {cce), bound-free 

photoionization (PI) and electron-impact ionization (El) cross sections {api,aEi), autoioniza¬ 
tion rate (AI) {Rai), and free-bound radiative recombination (RR) cross section {aRR), where 
api and apR are related through the Milne relation. 

The processes in Table 1 can be divided into two kinds, the inherent reactions inside the 
plasma and the reactions driven by an external radiation field. 

The inherent ones include the reactions between ions and electrons driven by collision, 
and spontaneous decay inside all ions. The solar corona plasma is believed to be dominated 
by these processes. In such low density plasma, the dominant processes are spontaneous decay 
and radiative recombination, whose rates are much higher than collisional decay and three-body 
recombination. As a result, the plasma departs from the local thermodynamic equilibrium, and 
can no longer be described by Saha and Boltzmann equations. The collisions between parti¬ 
cles cause energy exchange and state distribution changes. Usually charge exchange reactions 
between ions have negligibly low probabilities (Salzmann 1998). 

The rate per volume of electron impact excitation reaction is calculated by 
reactions 


cm^s 




( 2 ) 


where Ug is the electron density in plasma, and £c(cm^s ^) is the collisional excitation rate 
coefficient. The rate coefficient of collisional excitation is given by 

r+oo 


Sr.= 


r-\-oo 

/ f{v)acE{v)dv, 
J A.E 


( 3 ) 


where f{v) is velocity distribution of electrons, assumed to have Maxwellian distribution with 
electron temperature Tg, acE^v) is the collisional excitation cross section from Wj to Wy' at 
velocity v, and AE is the excitation energy. Thus, the CE rate coefficient expressed in term of 
incident electron energy is 

E, 


Sr = 


r+oo 

/ EacE{E)exp{- — )dE, 

JAE Ip 


( 4 ) 


irmpT^ 

where E is incident electron energy, rUe is the mass of electron, and (Jce{E) is the collisional 
excitation cross section calculated by FAC. 

The calculations of collisional ionization and radiative recombination rates have a similar 
form with CE. The Cl and RR rate coefficients (cnUs“^) are 

' E 


Xp = 


7rmgT| 


r+°° E 

/ Eaci{E)exp{- — )dE, 

JAE Ip 
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TZr = 


+ 00 

EaRR{E)exp{-—)dE, 
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and aRR^E) is the RR cross section 


TTTTleTg^ Jo 

where aci{E) is the Cl cross section from Xij to X*. 
from Xj+ij'/ to Xij. 

The inverse process to Cl is three-body recombination from to Xij and CE’s 

inverse process is collisional deexcitation from Xjj/ to Xij. Their rate coefficients are obtained 
by Detailed Balance Principle 

, 2 \ 3/2 


7^^ = - 


Vr = 
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-exp[-—)Xc, 
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9i,j" JJe 

It should be mentioned that the three body recombination rate coefficient, being dependent on 
two electrons, is proportional to 
reactions 


cm^s 






-JZt, 


( 9 ) 


and the unit of JZt is cm°s“ 

The spontaneous processes in plasma are radiative decay and autoionization. Their 
rates are directly given by FAC, and Rai{s~^). Radiative decay process is the 

mechanism of plasma emitting line spectrum. In the present work, we assume that the photons 
emitted by radiative decay are not reabsorbed. Autoionization occurs in case of doubly excited 
ions, and is possible only if the sum of the energies of the two electrons is higher than the 
binding energy of the ion. During autoionization, the energy, which is released by the inner 
excited electrons decays to a lower state, ionizes the outer one into the continuum. There are 
several ways to produce such highly excited ions, such as dielectron capture, photoexcitation 
and photoionization of inner-shell electrons, and the two inner-shell processes are especially 
important in conditions where strong helds exist. Dielectronic capture into doubly excited 
states is the inverse process of autoionization, and its rate coefficient (cm^s“^) is obtained by 
the detailed balance principle, 

, 2 \ 3/2 


Rd = - 
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The rate coefficient of dielectronic recombination is obtained when Eq.(lO) is multiplied by the 
branching ratio for radiative stabilization of the doubly excited state, Ai^/+ HRai) 
(Salzmann 1998). 

When the plasma is irradiated by a strong external radiation source, the radiative held 
will excite or ionize the ions causing the plasma gets into photoionizational collisional radiative 
equilibrium regime. In this case, the photoionization and photoexitation processes are not 
negligible and may dominate the charge state distribution. For example, the strong radiation 
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field emitted by an accreting compact object is believed to be the main ionization mechanism 
of the highly ionized low density gas aronnd it. 

The photoionization reaction rate per unit volume is calculated by 
reactions r+oo 


cm^s 


r+oo 

= NijRpi = Nij / np{hu)capi{hu)d{hu), 

JAE 


where Rpj is the photoionization rate(s“^), AE is the ionizing energy from Xij to Xj+ij/, hu 
is the energy of incident photon, apiijw) is the photoionization cross section, c is the speed of 
light, and np{hv)d{hv) is the density of photons having energy in the range [hv,hv + d{hvy\. 
For a black body radiation source having radiation temperature Tr, with energy intensity I{hu) 
(eV/(cm^-s- eV)) and dilution factor a, np{hu) = {al(hn))/(chn), then Rpj becomes 
f+°° I{hv) 

Iae 


Rpj — cr 


hv 


-api{hv)d{hv), 


( 11 ) 


The photoexcitation rate per volume from Xij to Xij/ is 

reactions at r t/t \ n 

= NijRpE = 


cnRs 


is the Einstein ^-coefficient The photo-excitation rate(s 


isotropical black body radiation source is 


= « nhX-^£AA. 




( 12 ) 

irradiated by a diluted 
(13) 


9j 2(fc!') 

Currently, RCF is assumes a blackbody radiator, thereby the photoionization and pho¬ 
toexcitation rate are 
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Rpj — Q. 


JaE p t. 
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QTr — 1 
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eTr — 1 
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2.3. Input Parameters 

These equations require four input parameters in RCF, which are radiation temperature 
Tr, dilution factor a, electron temperature Te, and electron density n^. Tr is the temperature 
of the blackbody radiation source, a stands for the attenuation of radiation between the source 
and the irradiated plasma, and it mainly depends on the opacity and distance. Te and n^ 
are the properties of the irradiated plasma, and they are sufficient to describe plasma under 
coronal equilibrium. However, when the external radiation held is important, all the four input 
parameters are needed. 

In experiments these parameters are measured directly or deduced indirectly from some 
measured values. However, depending on the experiment setup, some parameters cannot be 
obtained. For example, the experiment by Fujioka et al. (2009) with silicon target provided all 
four parameters with some uncertainties, whereas Tg does not have a dehnite value in the Sandia 
experiment by Foord et al. (2004) with iron target. In astrophysics, Tr is estimated by the 
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observed continuum spectrum, and a is roughly deduced from the distance between two celestial 
objects by the inverse square law. Usually, Tg and rig of the irradiated plasma are deduced from 
some characteristic spectral line ratios. For example, the ratios of resonance, intercombination 
and forbidden lines of He-like ions are important diagnostics of electron density and temperature 
(Porquet & Dubau 2000). The radiative recombination continuum (RRC) is also an important 
method to diagnose electron temperature of plasmas. 

3. Simulation of Sandia Photoionization Experiment 

Photoionized plasma is a special kind of NLTE plasma. It is widely observed in the 
universe, such as low density nebula near accreting X-ray source. Recently, this kind of plasmas 
were also produced in laboratory using high power laser (Fujioka et ah 2009) and Z-pinch (Foord 
et al. 2004). 

In this section, RCF is applied to the photoionization experiment at Sandia National 
Laboratory Z-facility (Foord et al. 2004). In this experiment, a 165 eV near-blackbody radi¬ 
ation source was created to produce a Ue = 2 ±0.7 x 10^®cm“^ plasma (Foord et al. 2004) in 
photoionizational collisional radiative equilibrium regime (Wang et al. 2011). A distribution of 
iron charge states was deduced from the absorption spectrum. A number of papers (Foord et al. 
2004, Foord et al. 2006, Wang et al. 2011, Han et al. 2013, Liang et al. 2014) tried reproducing 
the measured charge state distribution using different models and computer codes. All these 
works assumed a steady state photoionized plasma, which is also adopted by RCF. 

In this experiment, only two of the parameters needed in RCF are specihed, which are 
Tr and Ug. T^, is a disputed focus of the former works, and it spans from 70 eV to 150 eV in 
different models (Foord et al. 2006, Wang et al. 2011). However, a was not specified by some 
models (Foord et al. 2006), although it is an important parameter that controls the influence 
of radiative field on the plasma (Han et al. 2013). Fortunately, the experiment yielded an 
ionization parameter = 25 erg ■ cm ■ (Foord et al. 2004) at the peak of the radiation 
pulse. ^ is a parameter related to the radiation held ^ = IOtt^ J/ue, and for an isotropic blackbody 
radiation held ^ = IQnaT^/rie, where J is the mean intensity and a is the Stefan-Boltzmann 
constant. Therefore, the ratio between experimental value and theoretical value stands for the 
attenuation of the radiative held, ie. a = ^exp/^theo, which can be derived as a = 0.85% —1.76%. 

Figure 1 displays the charge state distribution of the iron photoionization experiment 
predicted by RCF and comparisons with the experiment values and some previous works (Foord 
et al. 2006, Wang et al. 2011). Seemingly, RCF produces the result closest to the experiment, 
and almost every ion is within the experiment uncertainties. The average charge state in present 
calculation is < Z >= 16.12, which agrees well with the measured value < Z >= 16.1 ±0.2. The 
input parameters used here are Tr=165 eV, Ue = 2 x lO^^cm”^, a = 1.4% and Tg = 150eU. 
Te=150 eV agrees with CLOUDY and FLYCHK, and some other works, such as NIMP (Rose 
et al. 2004) and Han et al. (2013). a = 1.4% is in the interval deduced above. 

6 



Fig. 1. The charge state distribution calculated by RCF (Solid) for Fe photoionization experi¬ 
ment (Scatters), and comparison with PhiCRE (Dash Dot Dot), GALAXY (Dash), CLOUDY (Dot), and 
FLYCHK (Dash Dot) (Foord et al. 2006, Wang et al. 2011). The temperatures are electron temperatures 
used by these models. 

A main reason of the differences among the codes in Figure 1 is the different sources 
of atomic data. GALAXY employs an average-of-con£guration approximation for electronic 
states, screened hydrogenic for both collisional and radiative processes, and Hartree-Dirac- 
Slater or Kramers cross sections for photoionization (Rose 1998, Foord et al. 2006). FLYCHK 
uses hydrogenic approximation to calculate energy levels and level populations (Chung et al. 
2003, Foord et ah 2006). Results of GALAXY and FLYCHK largely deviate from the measured 
one. The atomic databases of CLOUDY are accurate enough to be comparable with spectral 
emission data (Ferland et al. 1998, Foord et al. 2006), but there still are some obvious disparities 
between it and the experiment. The energy levels and spontaneous decay rates of PhiCRE are 
taken from the NIST database, and other rate coefficients are calculated by widely used formulas 
(Salzmann et al. 2011, Wang et al. 2011). 

The atomic data of RCF are calculated by FAC, which calculates all the atomic data 
by a fully relativistic approach based on Dirac equation (Gu 2008), and this single theoretical 
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framework ensures self-consistency between different parts. The configurations calculated with 
FAC of present work are listed Table 2, which include 1948 singly or doubly excited levels. For 
saving computation time, the maximum principle quantum number n is set to be 4. Because 
the states with K-shell vacancies have energies higher than 7 keV, which are much higher than 
the energies of photons and electrons under this experiment condition, K-shell is closed in 
FAC calculation. For ensuring the accuracy of present work, we compare the present FAC 
data with some literature values. Table 3 is the comparison of energy levels for 2s‘^2p^3l and 
2s2p^3l states of Fe^®^ between the NIST database (Kramida et al. 2014) and present work, 
and it shows excellent agreement (within 0.4%). Figure 2 shows the comparison of radiative 
decay rates (s“^) between present work and the available corresponding transitions on the 
NIST database (Kramida et al. 2014) of the four most abundant charge states. As shown, 
there are more than 80% of present data within 20% agreement with NIST database (Kramida 
et al. 2014). According to Eq.(14), the accuracy of radiative decay rates also guarantees the 
calculating of photoexcitation rates. For collisional excitation. Figure 3 is the comparison of 
cross section of Fe^^+ and Fe^®+ for transitions from their ground states to their Erst excited 
levels. The present FAC results agree well with those calculated with ICFT (intermediate¬ 
coupling frame transformation) R-matrix (Liang et al. 2009,Liang & Badnell 2010) and Dirac 
R-matrix (Aggarwal et al. 2003). Using these data, RCF successfully reproduces the experiment 
result, and we look forward to applying it to spectral analysis of laboratory or astrophysical 
plasmas in future works. 
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In Figure 4, we compare the influences of some important processes on the experimental 
results. Case A is the one presented in Figure 1, in which we included all the relevant processes 
in RCF and all the configurations in Table 2. 

Case B uses data without the doubly excited levels in Table 2. As can be seen, there 
are signihcant differences relative to Case A. 

In Case C both the autoionization and its inverse process, dielectronic capture, are turned 
off. Results of Case C are close to those of Case B, and both predict an average charge lower 
than that of Case A. This means that doubly excited states are non-ignorable in the current 
calculation of charge state distributions. In other words, autoionization is an important ionizing 
channel, and doubly excited states act like ladders to the next charge state in this experiment. 

In Cases D and E, we examine the influence of the radiation held on the Sandia experi¬ 
ment charge state distribution. First we note that there is a big difference between the binding 
energy of the two most important charge states in the experiment: the ionization energy of 
Fei5+ (Na-like) is 489.312 eV, and that of Fe^®+(Ne-like) is 1262.7 eV. As shown in Figure 5, 
a larger fraction of photons than electrons has sufficient energy to excite or ionize Fe^^^ and 
Fe"*"^®. In fact, we have found that the inhuence of electron collisional ionization is so small 


Table 2. The configurations used by Case A. 


Charge State Singly Excited 



Doubly Excited 


^gl4+ 

2s^2p^3s^ 

2s^2p^3s3p 

25^2/3/ 

2s'^2p^3s'^3p 

2s2p^3s^3p 


2s^2p^ZsM 

2s^2p^3p3d 

2s22/3/ 




2s^2p®3s4s 

2s^2p®3s4p 

2s^2p^3sAd 




2s^2p^3s4:f 





^gl5+ 

2s^2p^Zs 

2s^2p^3p 

2s^2p^3d 

2s^2p^3s^ 

2s'^2p^3s3p 


2s^2p^As 

2s^2p^Ap 

2s^2p^Ad 

2s22/3/ 

2s^2p^3s3d 


25^2/4/ 



2s^2p^3p3d 

2522/3/ 





2s2/3s2 

252/3s3p 





2s2/3/ 

2s2p^3s3d 





2s2/3p3d 

252/3/ 

^gl6+ 

2s^2p^ 

2s^2p^3s 

2s^2p^3p 

2s2/3s2 

2s2p^3s3p 


2s^2p^M 

2s2/3s 

2s2p^3p 

2s2/3p2 

2s2p^3s3d 


2s2p^M 

2s^2p^4:S 

2s'^2p^Ap 

2s2p^3p3d 

252/3/ 


2s^2p^Ad 

2522/4/ 




^gl7+ 

2s^2p^ 

2s2/ 

2s‘^2^3s 

2s2/3s2 

252/3s3p 


2s^2^3p 

2s^2‘^3d 

2s2/3s 

2s2/3p2 

252/3s3d 


2s2p^3p 

2s2/3d 

2s2244s 

2s2/3p3d 

252 / 3 / 


2s^2^4:p 

2s^2'^Ad 

2s^2‘^Af 




2s^2Hs 

2s^2^Ap 

2s^2Hd 




2s^2Hf 





^gl8+ 

2s'^2p^ 

2s2/ 

2/ 




2s^2p^3s 

2s2/3s 




^gl9+ 

2s23p3 






that its omission from the computations is hardly different from Case A. 

There are two radiation driven processes in the code, photoexcitation (PE) and pho¬ 
toionization (PI). PE is omitted in Case D, whereas PI is turned off in Case E. It can be seen 
from Figure 4 that the absence of PE signihcantly reduces the average charge state, he., pho- 
toexcited states play an important role in the charge state distribution. Case E is very close 
to Case A. The reason for this behavior is the large difference between the binding energies of 
and Fe’*'^®. On the other hand, there is a threshold between Fe^®+ and Fe^^’*', because 
neon-like Fe’*'^® has a closed shell stable configuration. 

Both of PE and PI have preference to ionize inner-shell electrons resulting in autoionizing 
states. However, according to the photon energy distribution in Figure 5 and the energy levels 
of main ions, the doubly excited levels seem to be more likely produced by PE process. Actually, 
when the PE channel to doubly excited states is shut down, the result is almost the same as 
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Table 3. Comparison of energy levels for 2s^2p^3l and 2s2p^3l states of Fe^®"^ between NIST and present work. 


Index 

Level 


NIST 

Present 

0 

2s^2p^ 

'^0 

0 

0 

1 

2s^2p^3s 

^P2 

725.2443 

723.9074 

2 


^Pl 

727.1388 

725.874 

3 

2s^2p^Ss 

^Po 

737.856 

736.5231 

4 


^Pl 

739.0537 

737.762 

5 

2s^2p^3p 


755.4915 

755.008 

6 

2s^2p^3p 

^D2 

758.9928 

757.7623 

7 



760.6095 

759.327 

8 



771.0614 

769.7896 

9 

2s^2p^ip 

^Pi 

761.7403 

760.5093 

10 

2s'^2p^3p 

^P2 

763.5529 

762.2596 

11 


^Po 

768.981 

767.8495 

12 


^Pl 

774.3073 

773.1159 

13 

2s^2p^ip 

^D2 

774.6855 

773.415 

14 

2s^2p^ip 


787.7224 

790.3189 

15 

2s^2p^M 

^Po 

801.4313 

800.4122 

16 


^Pl 

802.401 

801.3396 

17 


^P2 

804.211 

803.0749 

18 

2s^2p^M 


804.2644 

802.9189 

19 


^Fs 

805.0331 

803.6367 

20 


^F2 

817.5964 

816.2671 

21 

2s^2p^M 

^D2 

806.728 

805.3275 

22 

2s^2p^M 


807.8004 

806.4032 

23 



812.369 

811.2368 

24 


^D2 

818.4135 

817.0462 

25 

2s^2p^M 


818.9342 

817.4908 

26 

2s^2p^3d 

^Pl 

825.7 

825.4368 

27 

2s2/3s 


869.1 

867.1547 

28 

2s2p^3p 

'Pi 

892.55 

892.5898 

29 

2s2p®3p 

'Pi 

896.939 

895.3807 


Case D, which confirms that PE is the main pumping mechanism of doubly excited states. 
Case E indicates that PE+AI process wins the competition in ionizing Fe^®^, but for Fe^®^ the 
PI channel is important, too. 

In Case F, collisional excitation (CE) is omitted, and it gets a strange result. The 
deletion of CE does not lower < Z > as Case D, but rather the plasma is more ionized than 
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A 


NIST 


Fig. 2. Comparison of radiative decay rates between present work and corresponding available data on 
NIST from to Fe^^'*'. 

Case A. According to Figure 5, although the electrons have comparatively lower energy and 
cannot ionize the ions as effectively as photons, they still can excite the ions to singly excited 
levels by collision. However, the collision produced singly excited levels have smaller reaction 
cross sections with photons than the ground state and lower levels, namely, they are more 
difficult to be ionized by photons. Therefore, the electrons certainly would take part in the 
competition of reacting with the ground state of ions, and as a result reducing the ionizing 
efficiency of PI and PE+AI. What is more, when the CE channel from the ground states of 
ions to the singly excited states is shut down, it produces a result similar to Case F, which 
confirms the discussion above. So that, it makes sense that why < Z > rises when CE is shut 
down. 

In conclusion, RCF has a good agreement to the photoionization experiment results 
(Foord et al. 2004), and gives reasonable explanation for the charge state distribution. In 
the calculations of RCF, the charge state distribution of this experiment is a composite result 
of different atomic processes. The external field dominates the ionizations in the plasma by 
photoionization directly and photoexcitation plus autoionization indirectly. The transitions 
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Incident Energy (eV) 


Fig. 3. Comparison of collisional excitation cross sections. (a)Transition from 2p®3s ^S' 1/2 to 2p^3p ^Pi /2 
of Grey solid: ICFT R-matrix (Liang et al. 2009), black square: present FAC. (b)Transition from 

2s^2p® ^5'o to 2s^2p^3s ^P 2 of Fe^®+. Grey solid: (Liang & Badnell 2010), black square: present FAC, 
circle: Dirac R-matrix (Aggarwal et al. 2003). 

within any given single charge state can significantly affect the charge state distribution, and 
one of the interesting results of our computations is the role played by collisional excitation in 
this experiment, in which it reduces the total ionization rate by competing with PE and PI. 

4. Summary 

In this paper we introduced a new code, RCF, which is applied to plasma in NLTE 
condition, especially in the photoionization dominating regime. This code can calculate the 
level population, charge state distribution and spectra of a plasma in steady state. The atomic 
data source of this code is FAC, which is an easy to use and powerful software package to 
calculate various atomic data. The SEAC interface can provide all the atomic data needed in 
RCF, without any additional modifications. FAC is based on a fully relativistic theoretical 
framework, which ensures the accuracy and consistency of the atomic data. 

All the plasma processes and their inverse ones are related by the detailed balance prin- 
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charge state 


Fig. 4. The charge state distribution calculated by RCF with different data or without specific processes. 
Case A uses data with inner-shell holes (Solid), Case B use data without inner-shell holes (Dash).Case 
C (Dot), D (Dash Dot), E (Dash Dot Dot) use the same data as Case A, but there are no AI and DC 
processes in Case C, no PE process in Case D, no PI process in Case E, no CE process in Case F. The 
scatters are experiment values. 

ciple in RCF. As a result, in high density regime, the RCF generates results similar to the Saha 
equation with same atomic data. In other words, RCF converges to LTE approximation under 
the appropriate condition. In radiation dominant regimes, RCF gets a charge state distribu¬ 
tion which closely agrees to the results of the Fe photoionization experiment. A comparison 
is given to the results of other similar codes. We also discussed the influence of the various 
atomic processes to the charge state distribution of this experiment. Photoionization is not 
the only important ionizing channel, but the photoexcitation plus autoionization process are 
proved to be also signihcant. Although the electrons have comparatively lower energy than the 
photons, they still are important. The electrons can excite the ions to the levels which have 
small reaction cross sections for the photons, and the result is reducing the ionizing efficiency 
of the photons. 

The charge state and levels’ distributions are a prerequisite for the simulation of the 
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Fig. 5. The distribution of 165 eV blackbody photons (solid) and 150 eV Maxwellian electrons (dash). 

emission spectrum, and we showed that all the atomic processes may have signihcant effect 
under the appropriate conditions. In particular, in the analysis of X-ray spectrum from a 
compact object photoexcitation is an important pumping mechanism (Kallman et ah 2014). 
Porquet & Dubau (2004) also emphasized the influence of cascading decay from higher levels 
and collisional excitation on the line ratios in plasma diagnosis. According to the result shown, 
RCF is a reasonable code to get accurate distributions in steady NLTE plasma by including all 
the processes and using FAC data. We shall use it for spectrum analysis in astrophysics and 
laboratory of photoionizing and collisional NLTE plasmas in our further works. 
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